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Abstract 

We study a two-dimensional rotating Bose-Einstein condensate confined by an anharmonic 
trap in the framework of the Gross-Pitaevksii theory. We consider a rapid rotation regime 
close to the transition to a giant vortex state. It was proven in |CPRY3| that such a transition 
occurs when the angular velocity is of order £~‘^, with denoting the coefficient of the 
nonlinear term in the Gross-Pitaevskii functional and £ <C 1 (Thomas-Fermi regime). In this 
paper we identify a finite value Qc such that, if with Hq > He, the condensate is in 

the giant vortex phase. Under the same condition we prove a refined energy asymptotics and 
an estimate of the winding number of any Gross-Pitaevskii minimizer. 
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1 Introduction 

Since the first experimental realization of Bose-Einstein (BE) condensation in the 90’s, BE conden¬ 
sates and cold atoms in general have been extensively studied to investigate quantum properties 
on almost macroscopic scales. Among the typical features of BE condensates, one of the most 
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striking is certainly superfluidity, which has been studied in several experiments in the last years 
by putting the quantum system under rotation and observing its response (see, e.g., the reviews 
EillM]). Because of the quantum nature of BE condensates, the only possible change to the 
condensate profile due to the imposed rotation is the nucleation of isolated defects, i.e., quantum 
vortices. The generation of vortices has been observed in various experiments as well as the growth 
of their number when the angular velocity increases [MCWDl IRAVXKl ICHES) . For even larger 
angular velocities the number of vortices becomes so large that they hll the bulk of the system 
and arrange in a typical Abrikosov lattice |ARVK] . In presence of harmonic trapping the rotation 
can not be arbitrarily fast, otherwise the centrifugal forces would break down the trapping and 
system would eventually fly apart. On the opposite when the trapping contains some stronger 
confinement, e.g., some anharmonic potential growing faster that |rp for large |r|, regimes with 
much more rapid rotation can in principle be reached. Unfortunately so far a loss of coherence of 
the system has prevented the exploration of such regimes in the experiments [BSSDj . although a 
depression at the center of the trap has been observed for large angular velocities. 

However it has been predicted |CD[ ICDY^ IFe21 lEJS! IFB[ IKTUi iKBl IKFl IRT] that, besides the 
nucleation of vortices, other phase transitions should be observed in rapid rotating condensates in 
case of anharmonic confinement, with the occurrence of macroscopic defects or the transition to 
giant vortex states: when the rotational velocity gets very large, the centrifugal forces constrain 
the condensate in some thin annular region around a macroscopic hole and, if the rotation gets 
even more rapid, vortices disappear from the bulk of the system, which seems then to carry a huge 
circulation centered at the origin. 

Although BE condensates are many-body quantum system composed of a number of atoms 
ranging from few thousands to many millions, all the physical prediction about them are made by 
using an effective theory, the Gross-Pitaevskii (GP) theory, namely a one-particle approximation 
in which the energy of the system is given by a suitable nonlinear functional (see below). In 
spite of its simplicity the agreement with experimental observations is quite good, specially in 
the so called Thomas-Fermi regime, i.e., when the effective coupling becomes large. One of the 
major advantages of GP theory is the possibility of run very sophisticated and accurate numerical 
simulations [Danl IF.TSl IKTUj . See also the webpage http: //gpelab .math, cnrs. fr/, where one 
can find an efficient free code for simulations of the GP energy or dynamics developed by X. 
Antoine and R. Duboscq |AD11IAD^ . 

In the framework of the GP theory the energy of a two-dimensional 
physical units on the plane orthogonal to the rotational axis reads 

fp°hys['I'] = / dr {i |(V - zA,ot) + (U(r) - 

where Hrot is the angular velocity, Aj-ot := r = |r|, with r 

[—y,x)l\v\ is the unit vector in the transverse direction. The trapping 
of the form 

V{r) ■.= (1.2) 

with k > 0 and 

2 < s < oo, (1-3) 

i.e., the harmonic trapping is corrected by some anharmonic perturbation. Finally, we will focus 
on the study of the Thomas-Fermi (TF) regime e —>■ 0. The ground state energy of the system 
is thus obtained by minimizing the functional (HH) under the normalization constraint ||ik ||2 = 1, 
which amounts to require conservation of the particle number. Any minimizer is called condensate 
wave function and its modulus square, i.e., the associated probability distribution, is what can be 
observed experimentally. 

The range of validity of the GP description as well the derivation of the GP effective theory 
from the quantum mechanical description of a condensed Bose gas is an interesting topic on its own, 


rotating BE condensate in 

( 1 . 1 ) 
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which has been completely solved in the non-rotating case [LSYll ILSSY) . In presence of rotation 
on the other hand m contains a derivation of the GP functional, which is however restricted to 
bounded angular velocities and therefore not directly applicable to the case under discussion (see 
also [BCPY] for further results). However we will not investigate further such questions and take 
as a starting point the GP theory. 

The mathematical physics literature contains now a large number of works studying the behav¬ 
ior of the GP minimization problem in different asymptotic regimes of the angular velocity. If we 
restrict the discussion to trapping potentials of the form (EH), three phase transitions have been 
identified (see |CPRY3] for an extensive discussion or [GPRY4) for a more concise exposition), 
corresponding to three critical values of the rotational velocity. Here we briefly sum up the most 
relevant features of the physics of rotating condensates in anharmonic traps: 

• for small angular velocities Hrot, the rotation has no effect on the condensate wave function, 
i.e., the minimizer of coincides with the one in absence of rotation |A,TR,) : 

4 

• when the first critical speed Hd oc e^|log£:| is crossed, one observes the nucleation of 
quantum vortices, i.e., isolated zero of the condensate wave function |GR1) : 

a-2 

• if Hpot stays far from a second critical speed f 2 c 2 oc e , the number of vortices might 
increase but the profile of the condensate wave function is still close to the non-rotating 
one. Glose to 17^ it is possible to derive the explicit distribution of vortices [GRlj . which 
eventually cover the whole bulk of the condensate. In this regime one expects that they 
arrange in a regular (Abrikosov) lattice to minimize the interaction energy. This remains 
an open question although it has been proven that the vorticity is uniformly distributed 
|GPRY.^ : 

_ 5-2 

• for Orot oc £ 0+2 a first change of the macroscopic profile of the condensate is observed, due 
to the effect of the centrifugal forces. When the second critical speed Hcj is crossed this 
change has a dramatic effect since a macroscopic hole is created at the center of the trap. 
However the vorticity remains uniform in the bulk of the system |GPRY3] : 

• for very rapid rotations above flcj, the bulk of the condensate becomes essentially annular 
and its width shrinks as e —>■ 0. No further changes are however observed until a third critical 

4 ( 3 - 2 ) 

speed Hc 3 oc e is crossed. Then vortices are expelled from the bulk and the condensate 

behaves as if the whole vorticity was concentrated at the origin of the trap [GPRY31IGPRY5) . 
This is the giant vortex state that we plan to study in this paper. 

So far we have only discussed condensates in anharmonic traps of the type (na) but a lot of 
results are also available for other classes of trapping potentials. First of all the harmonic case has 
been extensively studied both in the physics and mathematical literature and, while there exists 
a first critical value of the angular velocity [IMll IIM2) corresponding to the occurrence of vortices 
and the behavior of the condensate for not too rapid rotation is similar to the one described above 
(vortex lattice, uniform distribution of vorticity, etc.), when the angular velocity approaches the 
harmonic frequency of the trap, some new physical features come into play and fractional quantum 
Hall states emerge |ABD[ lABNl ILSY2] . As we have already mentioned larger angular velocity 
are not allowed because the system would otherwise be no longer trapped. See however |Ka| for 
an alternative setting in which the trapping potential is suitably rescaled to reach fast rotation 
regimes. 

Even if we restrict to the anharmonic traps (ESI) there is an extreme case which is of certain 
interest, namely s = oo. Formally this corresponds to a confinement of the system to a two- 
dimensional disc of unit radius. Naturally one has then to provide suitable boundary conditions 
and both the Neumann [GDYll IGYl IGRY] and Dirichlet [GPRYl] cases have been deeply studied. 
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Indeed phase transitions analogous to the one described above has been found out even in this 
extreme case, although the nature of the third one is much more subtle. 

Let us now go back to the functional (HH) and introduce more convenient parameters: if we 
set 

Uphys := ~ (1-4) 

and obviously assume that flosc < I^rot, the trapping potential can be cast in the form 


V(r) = kr^ - (1.5) 

Since we are interested in exploring a regime in which both e —>■ 0 and Ophys —^ oo (or, equivalently, 
flrot oo), it is convenient to rescale units in the GP functional, in order to observe a non-trivial 
behavior [CPRY31 Sect. LA]: if one would naively minimize under the mass constraint 

Ijdtjl = 1, one would get trivially that the ground state energy diverges and the corresponding 
minimizer tends to 0 pointwise. The appropriate rescaling depends however on the asymptotics of 
I^phys and, in the regime we want to explore (very fast rotation), 


Ilphys > e “+G 

which leads to the rescaling (see |CPRY31 Sect. LA]) 

r = i?mx, ^'(r) = i?“V(x), flphys = An = fixe,?, 

where i?in stands for the unique minimum point of the potential (HH), i.e., explicitly 


:= 



( 1 . 6 ) 


(1.7) 


( 1 . 8 ) 


Under the scaling (113, the GP functional (EH becomes 

^pTy.[4^] = + (f - ^) , (1-9) 


with 

[ dx|il(V-zAo)V'l' + f^'W^(x)|V'P+e-'IV'lH 


( 1 . 10 ) 


The rescaled potential 

lU(x) := i (x« - 1) - 4 (x2 - 1) (1.11) 

is positive and has a unique minimum at x = 1, i.e., inf 3 ,gR+ IU(x) = IU(1) = 0. The rescaled 
angular velocity fl is related to the original physical quantities via 


n = ( 1 . 12 ) 

and condition dm) becomes 

(1.13) 

From now on we will focus on the analysis of the minimization of the functional (11.1011 on the 
domain 


:= inf 


We also set 


(1.15) 
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and denote by any minimizer, which is known to exist by standard arguments. In addition any 
' 0 ®^, which might be non-unique due to a breaking of the rotational symmetry and the occurrence 
of isolated vortices, solves the variational equation 

— ^ (V — iAn)^ 0^^ -I- + 2s~^ , (1-16) 

where the chemical potential (Lagrange multiplier) is fixed by imposing the normalization of 
0 GP; 

^GP^E^^+e-^f dx|0GP|". (1.17) 

As discussed in details in |CPRY3l Sect. LB], when fl the condensate has already 

crossed the second critical speed, i.e., its profile approaches a density supported on an annulus 
centered in the origin, whose inner and outer radii tend to 1 as e —?► 0. More precisely 10*^^ p is 
close in LP, p < oo, to the TF profile 

pTF(x) = i [pTF _ e^n^W{x)] ^ , (1.18) 

with /r'FF the chemical potential fixed by the L^-normalization of the function. A straightforward 
analysis shows indeed that is compactly supported and supp(p"'’^) = [xin,a;out] with |CPRY31 
Eq. (2.7)] 

a;out - a;i„ = C'(efl)“^'^^< 1, Xin/out = ^ + ), (1.19) 

as it can be proven by taking a Taylor expansion of W around a; = 1 in (ll.lSp and imposing the 
normalization. 

The vortex structure of 0*^^ is richer: being well above the first critical speed \ logej for 

the nucleation of vortices, the GP minimizer contains a very large number of vortices distributed 
all over its support. More precisely one can prove that the vorticity is uniformly distributed in the 
bulk of the condensate. As in |CPRY3l Eq. (1.42)], we denote by TT-buik C supp(p"^^) a suitable 
annulus {x j x< < a: < a;>} with = Xout/in + o((efl)“^/^). 

Theorem 1.1 f |CPRY^ Theorem 1.2]). 

If £~^ ^ ^ “C £~'^ as £ —>■ 0, there exists a finite family of disjoint balls {Bi} := {R(xi, p0} C TT-buik, 
i = 1 ,..., N, such that 

1 . g.<0 ( 12 - 1 / 2 ), ^^2 < (1 + (ef2)2/3)-i; 

2 . > 0 on dBi, i = 1 ,..., N. 

Moreover, setting dt := deg{0‘^i^, dBi} and defining the vorticity measure as p := ‘2TrEf=i di6(x — 
X,;), then, for any set S C TT-buik such that j55j = 0 and j5j ^ fl-ij log(£:^f2)j^ as £ ^ 0, 


l^jS) 


( 1 . 20 ) 


The inner region {x € j a: < a:in} is presumably also filled with vortices but, because of the 
exponential smallness of 0 ^^ there, the vortex structure in that region is practically inaccessible. 
An important condition contained in Theorem 1 1.1 1 is the request 


n<^£-\ 


The reason is that at angular velocities of that order the proof of Theorem 11.11 might fail due to 
the occurrence of a further phase transition, i.e., the transition to a giant vortex state. This paper 
is precisely devoted to the investigation of such a transition. 
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From the heuristic point of view it is quite simple to explain why one should expect a change in 
the vortex structure when fl ~ £~'^: from energy considerations it is easy to see that the average 
size of the vortex core, i.e., the radius of the region around a vortex when is substantially 

far from , is of order . The width of supp(p'^^) is on the other hand of order 

and the two quantities are clearly of the same order when ~ Hence it must happen that for 
O = with Hq a large enough constant, the vortex core becomes larger that the bulk of the 

condensate, i.e., vortices can no longer be accommodated in supp(p"^^). A non trivial phase factor 
of '0^^ is however needed in order to compensate the effect of the rotation but, because no vortex 
can occur in the bulk of the condensate, all the vorticity should get concentrated in the inner region 
where 0^^ is exponentially small. In fact when this occurs it is impossible to distinguish from the 
energetic point of view such a state with vortices distributed in the inner hole from a giant vortex 
state of the form n G Z. 

Notice that although this might seem to suggest that the rotational symmetry is restored, 
such a phenomenon never occurs as proven in |CPRY3l Theorem 1.6]. However the GP energy 
is expected to be well approximated above the critical speed for the transition to a giant vortex 
state by a one-dimensional energy functional obtained by evaluating on functions of the form 
/(x)e™’^. In fact by some very simple observations one can show that n = [HJ(1 -|- o(l)), where 
[ • J stands for the integer part. Let us now fix the angular velocity to be 

0 =^, ( 1 . 21 ) 

with Hq a positive constant. Concerning the giant vortex regime, the main results proven in 
[CPRY8) are stated below. We denote by Abuik a suitable annular layer around x = 1 containing 
the bulk of the condensate (see next (12.81) for a precise definition). 

Theorem 1.2 ( [CPR.Y31 Theorem 1.3]). 

If El is given by (ll.2ip . there exists a finite constant CIq such that for any EIq> EIq, no minimizer 
0®^ has a zero inside Abuik if £ is suffciently small. 

Theorem 1.3 f [CPRY31 Theorem 1.4]). 

If H is given by (11.211) with Hg > Hg as in Theorem 11.ill then as e (fl 


ijGP 


min 

Il/ll2 = l 




+ 0(llogel9/2). 


( 1 . 22 ) 


The first result, although being a consequence of the energy asymptotics p.22l) . is the most relevant 
one, since it shows the occurrence of the giant vortex transition for angular velocities of order e“^. 
The precise mathematical statement is a pointwise estimate in the bulk region of j0®^ j in terms of 
a strictly positive function, i.e., the minimizer of the functional appearing on the r.h.s. of (11.221) : 
since the latter is bounded from below by a positive constant in the bulk and the difference is 
pointwise small in e, also 0*^^ can not vanish there. 

For the analysis of the present paper it is very important to remark that both results hold true 
if the angular velocity is expressed by (11.211) with fig large enough, namely no precise estimate is 
derived there on the sharp value for the transition (see also Remark |2.31) . We indeed expect that 
the giant vortex structure appears as soon as El becomes (asymptotically) larger than 

(1-23) 

for some explicit value He. In this paper we will indeed investigate such a question and exhibit a 
hnite value He which is a good candidate for the sharp constant. Actually we are going to see that 
such a constant is a solution of some algebraic equation (see (12.710 involving quantities relative to 


^We use here polar coordinates [x, tt) G IR+ X [0, tt) on the plane. 
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a limit problem independent of e. Although we have not proven it yet, we do expect that next 
(12.71) has a unique solution, thus providing the sharp value of the critical velocitjU. 

We outline here the structure of the paper. Next Section contains the main results, i.e., the 
identification of the explicit value of the angular velocity for the transition to the giant vortex 
state, together with an asymptotic expansion of the GP ground state energy which is actually 
on the main ingredients of the proof of the above mentioned result. We also show that, as in 
|CPRY31 Theorem 1.5], one can deduce a (better) estimate of the total winding number of any GP 
minimizer. 

Sections [3] contains some preliminary estimates and a detailed analysis of the effective 
functionals that will play a significant role throughout the proofs. In Section 13.41 we prove 
the main properties of the cost function and in particular its positivity, which is the main 
mathematical tool used in the proof of the giant vortex transition as in several other works 

|( IPHYI 1ICPW^ l( :PHY3l [CFn ICW] . 

Sections m and [S] are devoted to the proofs of the main results: we first (Sectional) obtain the 
asymptotic expansion of the GP energy by comparing suitable upper and lower bounds and then 
(Section [S]) use such a result to deduce the pointwise estimate of showing the absence of 

vortices in the bulk. 

Notation: In the asymptotic analysis e —>■ 0 we will often use the Landau symbols: given a 
positive function g, we say that 

• / = 0{g) (resp. = o{g)), if lime^o \f\/g < C < oo (resp. = 0); 

• f (X g, whenever lime_>.o \ f\/g = C, with 0 < C < oo; 

• if / ^ 0, / -C g is synonimous of / = o{g) and f ^ g simply means that g ^ f ■ 

Sometimes we will use the notation 0{ \ log£|°“) to indicate a quantity of order | log£|“ for some 
finite but possibly large a. Since such a quantity will typically appear multiplied by powers of £, 
the explicit value of a will be irrelevant. 

We denote by Bg{x.) any two-dimensional ball centered in x and with radius g and by [xj the 
integer part of the real number x. The symbol C will stand for a finite constant independent of £, 
whose value might change from line to line. 

Acknowledgements: The authors acknowledge the support of MIUR through the FIR grant 
2013 “Gondensed Matter in Mathematical Physics (GOND-MATH)” (code RBFRI3WAET). 

2 Main Results 

The first non trivial observation to improve the results proven in |GPRY3| is that instead of making 
a special choice of the giant vortex winding number ([IIJ in [CPRY3] 1. one might try and optimize 
w.r.t. such a parameter, so obtaining a better candidate for the giant vortex state. This leads to 
consider the functional obtained evaluating on a giant vortex ansatz f{x)e™^ and minimize 
w.r.t. both / and n to find out the optimal giant vortex phase, i.e., explicitly 

£f[g\ = J dy (I -b £^ 2 /) \yg{y)\^ + Ufi{y)g^{y) + e^y^v{y)g^{y) + ^g^{y)'^ , (2.1) 

where we have set for convenience n = + /3 and exploited the exponential fall off of '0®^ to cut 

the tails \y\> g x \ log£|. The spatial coordinate has also been rescaled around |x| = I by setting 

^Strictly speaking in order to show that Qc is the sharp value for the transition one should also prove that, below 
Qc, vortices are still present in the bulk of the condensate, as done in m for hard anharmonic traps. We will come 
back to this question later. 
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a; = 1 + e^y. The potentials Ujs and v are obtained via a Taylor expansion of W(x) around a: = 1 
and to the leading order in e are simply given by a shifted quadratic potential (see (13.81) and (13.911 
for their explicit expressions). We remark however that in C/g, the parameter /3 always appears 
multiplied (at least) by e^, so showing that the correction is only lower order. 

Setting := inf||j^||^i and denoting by the corresponding minimizer, which can be 

proven to exist and be unique (up to multiplication by a phase factor) (see Proposition 13.111 . one 
can subsequently minimize w.r.t. /3 G K, obtaining the energy , an optimal phase /3* and a 
density g*, i.e., 

■■= mm£;f = E^ = (2.2) 

In Subsection 13.31 we will prove that = 0(1), so that, by the above argument, one expects 
the functional to be close in the limit e —>■ 0 to the following simplified giant vortex functional 


£^''[9]= f dy{^{g'f + ^y^g^ + ^g'^}, 

with ground state energy E^'' and minimizer gg^, i.e.. 


E^'':= M £^[g]=£^[gg^i 


where 

:= {g e H\R) \ yg e = ij . 

Here we have denoted for short 

OL := Hoa/s 2. 

The minimizer gg^ solves the variational equation 


- ig" + da-^y^g + £g^ = 


(2.3) 


(2.4) 


(2.5) 

( 2 . 6 ) 


where = Es^ + ^ hsAr 

We are now in position to introduce the explicit value of the constant £1^ appearing in the 
critical value of the angular velocity Hcg, which can be expressed in terms of the critical quantities 
associated with the effective one-dimensional functional f®'' and, specifically, ggv and /x®'': we 
denote by He the largest solution of the equation 


Ho — 


4 

s 2 





(2.7) 


where the r.h.s. depends on Hq through y^'' and gg-y. The existence of such a solution is proven in 
Proposition 13.Ill Note that thanks to the estimate ||5gv||^ < tt^®'' (see (I3.18p l. He > 0. 

Before stating the main result of this paper, we have to define more precisely the region we 
identify with the bulk of the condensate: we set for any a > 0 


^buik := {x G I ggv (^5^) > |loge| “} , (2.8) 

and observe that by the exponential decay proven in Proposition 13.31 ik) ~ ^ 

i.e., it certainly contains the bulk of the system. 

Theorem 2.1 (Absence of vortices in Abuik)- 

//H = Ho/e"^ with Hg ^ ^ ^ Oj ih/GTi TLO (jP TTlZTlZTTtZZGT' ijj COTVtdVTiS 'VOT't'lCGS %Th 

preciselyfor any x G ^Ibuik 

|x^GP(x)| = (i + 0(ei/2|ioge|-)). 


(2.9) 
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Remark 2.1 (Giant vortex structure). 

The pointwise estimate (|2.9I1 suggests that is approximately radial within .Abuik- As already 

mentioned, this does not imply that the rotational symmetry is restored, since one expects that 
is far from being radial in the inner region x < Xin, where several vortices should presumably 
be distributed more or less uniformly. In any case no GP minimizer is invariant under rotations 
if n is large enough \CPB,YIA Theorem 1.6], e.g., in the giant vortex regime. 

Remark 2.2 (Third critical velocity). 

In Provosition \,‘l. 1 1\ we will prove that the equation (Ell) has a solution. Although we do not prove 
it, we strongly believe that such a solution is in fact unique and identifies the sharp constant in the 
value of the third critical speed. 

More precisely Theorem \2.1\ indicates that above Oc/e^ the system undergoes the phase transition 
to the giant vortex state and the bulk of the condensate becomes vortex free. Hence 

( 2 . 10 ) 

We actually expect that flea = which obviously requires to prove that the solution to (E3 

is unique. In addition one should also prove that for slower rotations vortices are still present in 
the bulk of the system. We plan to attack such a problem in a future work, but here we want to 
stress that the negativity of the cost function (see next Section HE) for rio < lie Is a very strong 
indication that vortices are indeed convenient in this case and thus the sharp value of the critical 
speed is precisely 

Remark 2.3 (Comparison with |CPRY3) 1. 

We want here to discuss in more details the comparison between Theorem \2.1\ and the analogous 
result proven in ICPB.YIA Theorem 1.3]: in principle, one could indeed derive an estimate of the 
threshold Hq for the transition to the giant vortex state there and then it would be natural to compare 
it with the explicit value found here. However we provide here some heuristic arguments showing 
that such a comparison is actually not needed (see however next Bemark \3.1\ for futher details). 
First of all an explicit estimate of Hg is not an easy task to achieve, due to the proof structure in 
\CPB.Y3] : the result proven there is indeed obtained through an asyptotic analysis as fig —t oo and 
one should then estimate all the coefficients of the error terms appearing in the formulae. Such 
quantities ultimately depends on the pointwise estimate of the difference between the giant vortex 
profile and the ground state of the harmonic oscillator given in \CPBY1A Proposition 3.5], which 
is not explicit at all. 

However, even assuming that one could obtain a sharp value Hg, there are strong reasons to believe 
that, unlike flc (see also the previous B,emark \2.2\) . it can not be the coefficient of the critical speed. 
First of all the condition ilg > fig guarantees the positivity of the vortex energy cost in \CPB.Y3f 
fBemark \2.2\) and therefore CIq > ^c- Moreover, as explained in ]CPRY3] (see also \CPBY2}] ]. 
when fig —^ oo, another transition takes place, i.e., the condensate density profile goes from a 
TF-like shape (11.181) to a gaussian function minimizing some suitable harmonic energy. The key 
fact is that such a transition is expected to take place after the giant vortex one. Indeed here we 
show that, for finite fig, when the profile change has not yet occurred, the condensate is already 
in a giant vortex state. On the opposite, a quiek inspeetion to the proof in \CPBY3] reveals that 
the transition to the giant vortex is proven there by imposing that the profile is already gaussian. 
Hence any so obtained threshold value can not be meaningful. 

Remark 2.4 (Giant vortex density). 

We have formulated the pointwise estimate (1^ with pgy, but an analogous statement holds true 
with pgv replaced with g*. The error in (12.911 is indeed so large that one can not appreciate the 
difference between the two reference profiles (see Provosition \3. .91) . Let us stress however that the 
use of 5 * as a reference profile in the proof is on the opposite crucial to obtain the result (compare, 
e.g., the asymptotics (12.111) and (I2.12 |) ). 
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The absence of vortices proven in Theorem 12.11 and the pointwise estimate of follows from a 
refined result about the energy asymptotics in the same regime, that we state in the following 

Theorem 2.2 (Energy asymptotics). 

// n = with flo > lie as £ —t 0, then 


pgv 

= ^ + 0 ( 1 ). 


( 2 . 11 ) 


Remark 2.5 (Energy expansion). 

The leading term Ef' je'^ contains the main energy contribution due to the inhomogeneity of the GP 
profile together with the subleading kinetic energy of The absence of vortices in .Abuik can 

be read in the very small remainder term 0(1). It is indeed interesting to compare (12.111) with the 
analogous result \CPRY3l Theorem 1-4], where the error term is much larger, i.e., 0(| loge|®/^), 
in addition to the fact that the result proven there holds true only for Qq large enough. 

Notice however that the coefficient of the leading term Ef^ still depends on e, through the 
boundaries of the integration domain as well as the optimal phase /3* and the potential If one 
wanted to extract a proper asymptotic expansion then the natural statement would be 

EGP = :^ + 0(|loge|7), (2.12) 

with a much worse error term. 

Thanks to the pointwise statement (ES), one can deduce that does not vanish on .Abuik- In 
particular for any R = 1 + > 0 on DBr. Hence it is possible to define the winding 

number of on OBr for any such R. A consequence of the energy asymptotics and the estimate 
(12.9|) is thus the following 

Theorem 2.3 (Winding number). 

Let H = r2o£“^ with IIq > 11^ and R be any radius such that i? = 1 + 0{£^) as £ —^ 0, then 


deg = ^ + C’(l). 


(2.13) 


Note that the combination of the above result with the proof of the rotational symmetry breaking 
given in |CPRY3l Theorem 1.6] implies the presence of vortices in the inner hole region where 
is exponentially small. 


2.1 Heuristics 

Before discussing the proofs of the main results, we briefly expose the proof strategy from a heuristic 
point of view, i.e., not tracking down the error terms and neglecting most technical points. As 
usual the main result about the behavior of the condensate wave function is deduced from the 
energy asymptotics (12.1111 . We thus focus on such a proof. 

Most of the relevant features of a fast rotating Bose-Einstein condensate were already discussed 
in details in |CPRY3) and recalled in the Introduction. Here we take as a starting point the effective 
functional m which is expected to provide the leading order term in the energy asymptotics in 
units £“^. Note that the ground state energy of always provides an upper bound to E^^ for 
any integer phase, i.e., whenever H + /? € Z. Actually the same upper bound can be proven to hold 
true up to some small error term even if H + /? is not an integer (see Section ST]). Hence we can 
neglect the upper bound part of the proof and discuss only the lower estimate to E^^. 

A preliminary step which is already described in details in |CPRY3) is the restriction of the 
integration in to the bulk of the condensate, i.e., to an annulus centered in the origin with 
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radius ~ 1 and width O(e^). This can be done by exploiting the exponential decay of "0®^ 
outside. From now on we will then assume that the integration in x is restricted to the annulus 
|1 — a;| < 0 {e^\ log£|). 

The main steps in the energy lower bound are then the following: 


1. optimal giant vortex phase and profile: we minimize w.r.t. to G K and obtain a 
minimizing /3* and an associated density g*. It is crucial to observe that such a minimization 
yields an additional equation involving g*, which is in fact nothing but the vanishing of the 
first derivative of Ef' w.r.t. (3. Such an equation will play a crucial role at point 4 below; 


2 . 


splitting of the energy: using a technique introduced in [LM] , which is now rather standard, 
we decouple ■0^^ = exploiting the variational equation satisfied by 

5 *, we obtain 


GP E^ , £[u] 

£4 ^27r£2’ 

with u essentially minimizing the reduced energy functional 


(2.14) 


£[u\ = J (bcgl |i|Vu|Va(x) • j„(x) + , (2.15) 


where the “magnetic potential” a depends on 17 and /I* and is the superconducting current 


j„(x) = §(uVm*-u*Vm). (2.16) 

Completing the lower bound means to show that £[u\ is positive; 

3. hydrodynamic estimate: we note that the “magnetic potential” is divergence free and there¬ 
fore it exists a potential function F{x) such that 2gJ(a:)a(x) = —V^E{x). This trick was 
first used in [CBY] in the context of the GP theory for rotating condensates. For later appli¬ 
cations to the GL function see also [CEllEEs]. We can thus integrate by parts the second 
term in (12.151) obtaining 

y'dxF(x) curl(ju). (2.17) 

At this stage we observe that since /3* = 0{1) and it appears in (12.11) always multiplied by 
£^, a good approximation of the functional £p^ can be obtained by taking the limit £ —>■ 0, 
which yields the functional (12.31) . with ground state energy and minimizer gg^- We can 
also replace F{x) with its limiting counterpart F^'^[x), which is in fact a negative function. 
The last step to estimate (12.171) is to use the trivial inequality |curl(j„) | < |Vit| and the 
negativity of F^'' to get the lower bound 

£[u] > y dx {yi, + FS-) \Vu\\ (2.18) 

where we have also dropped the last positive term in (12.151) : 

4. positivity of the cost function: the above lower bound suggests that any topological defect 
of u should carry an energy cost given by the cost function 

K^'^ = \gl. + F^. (2.19) 

Positivity of such a function in the bulk would then imply that vortices are not energetically 
favorable anywhere in the condensate. This is turn can be proven by direct inspection of the 
function itself. First we observe that both g^^ and F®'" are radial functions and we therefore 
change coordinates x = 1 -I- e^y, so that in the new variable y the bulk of the condensate is 
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basically the whole real line. In the new variable the explicit expression of (that we still 
denote by F^'') is 

poo 

F^^y) = -2no dttgl^it). ( 2 . 20 ) 

Jy 

Notice that by symmetrjH of g^v, F^''{—oo) = 0 and F^''{y) < 0 for any y € R. The cost 
function is therefore smooth and i^®'"(±oo) = 0, so that, if it becomes negative, it must 
have a minimum. The derivative of it'®'' can be easily computed 

= 9gv{y)9gv(.y) + ‘^^oyylviv), (2-2i) 

so that, by strict positivity of ggv, at any critical point yo for isT®'', one has 


y'gviyo) = -2fioyo3gv(2/o)- ( 2 . 22 ) 

Now using the variational equation for g^^ and manipulating the expression (12.201) of the 
potential function, it is possible to show that the cost function can be equivalently rewritten 
as 


X®'' = 


^oy 


a 

ira 


0 _ 2 
2 




5gv 


fio / 2 

o" 9gy 


and, inserting the condition (|2.22p satisfied at any minimum point tjq of , we get 


K^^'iyo) = 


1 “1“ 1) 2 
2 7+ 


2 I “^0 2 / \ 

^0 2 dgyiVo) 

7ra^ ^ 


5 gv( 2 /o)- 


(2.23) 


(2.24) 


Using the parity of ygv as well as the variational equation, one can prove that the quantity 
between brackets on the r.h.s. of the expression above is positive if and only if it is positive 
at the origin (see Proposition 13.101) 










> 0 


Uo > Uc. (2.25) 


Once the energy asymptotics is proven, the pointwise estimate of which allows to exclude 

the presence of vortices in the bulk for Uq > Uc, is a simple consequence: putting back the positive 
term we have dropped in the lower bound, one first obtains an estimate of the region where |u| can 
differ from 1. Then combining this with an estimate of the gradient of u, one gets the result. 

It is worth mentioning at this stage a technical difference with previous approaches. Indeed in 
|CPRY3) two potential functions were actually used instead of one, in order to get rid of boundary 
terms coming from the integration by parts described at step 3 (see the discussion in |CPRY31 
Sect. C]). Here on the opposite we are able to use only one potential function by estimating in a 
more refined way the boundary terms (compare, e.g., with next (14.161) 1. As in [CPRY3) we also 
exploit the symmetry properties of the profile y*, which is to a very good approximation invariant 
under reflections w.r.t. the origin. 


3 Preliminary Estimates 

Here we collect some useful technical results as well as the main properties of the effective function¬ 
als involved in the analysis. An important piece of information is contained in Section 13.41 where 
we prove the positivity of the cost function. 

^Unlike ggv, the profile g* is not exactly symmetric, but F satisfies analogous properties thanks to the optimality 
condition of /3*, i.e., the additional equation involving g* and which was mentioned at point 1. 
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3.1 Giant Vortex Functionals 

We start by describing the derivation of the functional ( 12 . 11 ) from the GP energy. As anticipated in 
Section[2]the idea is to evaluate the energy of a trial state of the form /(a;)e*"’^ in polar coordinates 
X = (x, i9) and with n = + /3. In addition we assume that / is real as it will be for any giant 

vortex profile. The result of a rather simple computation is 


poo 

£:GP[/e™’’] = 27r / dx x 
^0 


i|V/|" + l !2 



1 

n j 


W{x) 


P + 



(3.1) 


Exploiting the exponential smallness of outside of the bulk of the condensate proven in 
|CPRY3) and recalled in next Pror)osition l3.31 we can restrict the integration domain to the annulus 


Ar, := {x e : |1 - x| < e'^r]} , ry := loge|, 


(3.2) 


where 770 > 0 is an arbitrary finite constant and the prefactor in the definition of 77 has been chosen 
of that form for further convenience. Thanks to (13.191) 

{x) = 0{£°°), for any X ^ A,,, (3.3) 


and the restriction is thus well motivated. In addition we will also see that a similar estimate 
holds true for any giant vortex profile. In terms of the one-dimensional functional (EH) we are 
then integrating in the interval [1 — £^ 77,1 + e'^rj] and a change of variable is called for: setting 

x = l + g{y) = VTre /(I -h e^y) (3.4) 

so that g is normalized iifl ■= L^([— 77 , 77 ], (1 -|- e^y)dy), we obtain the energy 


£''''{ 9 ] = -, dy{l + s^y)U{g'f + 

^ J —77 ^ 

i (1 -k £^y - ) +W{l + £^y) 

Now we expand IT(1 -I- e^y) in Taylor series around y = 0 to get 
W{1 + e^y) = 


g^ + ^g^\. (3.5) 


(3.6) 


where if{y) = 0{y'^). Using this fact we can rewrite the potential in (13.51) as (recall that o? = 

+ 2 )) 






1 (\+e-^y) 

with V independent of (3 and of lower order w.r.t. io Up. Explicitly 

1 /..2 


= Up{y) + £ y y(.y), 


Uy{y) ■■= 

v{y) := 


“ -y^ - 2Qo£^I3y - ^oe'^Py'^ + 


(l + e2y)"V2" " " 2 

U§(s + (s-l)e2y) ^ (s-l)(5-2)fl§ ^2^2 

t;—-^- R -^—5-^5 (y) • 

(l-ke^y) 6 


Some trivial estimate using the Taylor expansion EH) implies that for y G [~y, y] 

Up{y) = \a^y'^ + O (e^(l -k |/3|)y -k , 


(3.7) 

(3.8) 

(3.9) 

(3.10) 


^We set in fact := LP([—r;, 77], (1 + e^y)dy) for any 1 < p < 00. 
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which shows that, if, e.g., /3 is uniformly bounded in e, the potential Up[y) is harmonic up to 
corrections of higher order in e. Alternatively one can think of Up as a shifted harmonic oscillator 
by writing 

Up{y) = ia2 + 0{e^\P\y^+e^f3^). (3.11) 

In fact, since the optimal value of /3 we are going to choose is 0(1), both representations are 
equivalent since the shift will be O(e^). Concerning the rest v{y) one trivially has the upper bound 

Hy)\<Cn,+0{e^r,), (3.12) 

for y G and with a finite constant Cho- The rest in the above expression is a consequence 

of the bound \‘f{y)\ < C\v\^^ which follows from the Taylor expansion (13.61) . 

In conclusion we have recovered the expression m, i-e., 

£f[9\ = J dj/ (1 + er^y) ( 5 ')^ + Up{y)g^ + e^y^v{y)g^ + . 

We now discuss the minimization of such a function w.r.t. g and for that purpose we have to 
identify the proper minimization domain, i.e., 

:= e H^{-'n,r]) \ g = g*, \\g\\i^2 = l|. (3.13) 

The ground state energy of £|'' is defined as 


Notice that the assumption g = g*, i.e., reality of the argument, does not imply any loss of 
generality because the ground state can always be chosen real (see next Proposition). 

Proposition 3.1 (Minimization of £f')- 

There exists a minimizer gp £ of (12.11) that is unique up to a sign, radial and can be chosen 
strictly positive. In addition gp £ C°°{—rj,g) and it solves the variational equation 

- hp - 2(i+e^y) 9'p + Ufii.v)9fi + £^v^v{v)9fi + I9I = Tp9y (3.15) 

with Neumann boundary conditions g'p^ztr]) = 0 and yp = + ^|| 5 / 3 ||^ 4 - 

Finally gp has a unique maximum point at yp and it decreases monotonically anywhere else. 

Proof. Existence and uniqueness of the minimizer follow from strict convexity of the functional 
with respect to the density p = g^. The variational equation ()3.15ll is satisfied at least 
in weak sense. Then one deduces the strict positivity of gp noticing that it is actually a ground 
state of a suitable one-dimensional Schrddinger operator. The equality for p,p follows integrating 
the (I3.15|) and recalling the fact that gp has T^-norm equal to one. Finally a trivial bootstrap 
argument allows to deduce smoothness of gp and therefore that (I3.15|) is solved in a classical sense. 

The only non-trivial result is the one about the existence of the a single maximum point for 
gp. However it follows from the property of the potential Up{y) + e^y^v{y): going back to the 
expression of the potential in (13.11) . one can easily compute, with a; = 1 -|- e^y, 

9 [Upjy) +£'^y^v{y)] _ J_ 

dx 



which vanishes at a single point j/pot, i-e., where 


1 + e^2/pot - (1 + Tif) 


= 1-b 
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The Taylor expansion also shows that 

~ (s+2)no 

The monotonicity property of (7^ can then be obtained by a simple rearrangement argument (see, 
e.g., [CPRYll Proposition 2.2]): since the potential has a single maximum point, if gp had more 
than one maximum besides yp, one could move mass from the further maximum to the minimum 
in between and lower the energy. Since gp is a minimizer one gets a contradiction. □ 

Another effective one-dimensional functional which is going to play an important role in the analysis 
is (12.3L i.e., the formal limit e —>■ 0 of , assuming that /? = o(e“^): 

^^'^[9]= f dy{lig'f + ^y^g^ + ^g^}. 

J M 

The minimization domain is in this case given by 

:= {g GH\R)\g = g*, = l} , (3.16) 

and the ground state energy will be denoted by = infgg®gv £^''[g\. 

Proposition 3.2 (Minimization of £^). 

There exists a minimizer g^^ € of (E31) that is unique up to a sign, radial and can be choose 
strictly positive. In addition g^v € G°°(K.) and it solves the variational equation 

- hgv + ^y^9sv + (3.17) 

with 7x8'' = £’8'' -h ^llffgvllt 

Finally g^v is even w.r.t. the origin and has only one maximum at y = 0, which fulfills the 
inequality 

ffgv(O) = IISgvIlL < (3-18) 

Proof. See the proof of Proposition 13.11 Parity of g^, is a trivial consequence of the parity of the 
potential. The inequality (13.181) follows from direct inspection of the variational equation (13.171) : 
at any maximum point g",, < 0, which immediately implies the result. □ 


3.2 Estimates of the Gross-Pitaevskii and Giant Vortex Profiles 


In this Section we collect several technical estimates of the profiles involved in the discussion. Such 
estimates will play a key role in the proofs but can be typically obtained by standard techniques 
in functional analysis. 

We start by recalling a result which was in fact proven in |CPRY31 Propositions 3.1 and 3.2]: 
let rjo be the parameter appearing in the definition (13.21) of Ari, then 


Proposition 3.3 (Exponential decay of 

If n = £lo/e^, there exists two finite constants c,C > 0 (independent of go) such that, for any 

X c) , 




max 


e 4 , exp 




(3.19) 


In particular the above result implies that by taking go large enough we can make arbitrarily 
small outside Ar/. This fact will be crucial in restricting the computation of the GP energy within 
A,j. Notice also that as soon as |1 — a;| ^ e‘^\ loge|, = 0{e°°). 

Let us now focus on the giant vortex profiles. Before stating the main technical estimates we 
first formulate a simple preliminary bound on the giant vortex energy : 
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Proposition 3.4 (Preliminary bound on 
If 13 = O as e —>■ 0, then 

£;|" = G(1), tip =0(1). (3.20) 

Proof. Since is positive (compare with (I3.5|) l. it suffices to prove a suitable upper bound: to 
that purpose one can simply evaluate the functional on the ground state of the one-dimensional 
harmonic oscillator with frequency a. The result on tip follows from the trivial estimates < 
tip < 2Ef. □ 

The giant vortex profile gp decays exponentially for large \y\ and one can actually show that this 
decay captures the correct asymptotics of gp: 

Proposition 3.5 (Pointwise estimates of gp). 

If (3 = O (e“^) as e: —>■ 0, then there exists a finite eonstant C such that 

9 p{y) < , so that gp{±g) = O (e’'“). (3.21) 

If (3 = 0(1) then there exist two finite constants Ci,C 2 > 0 such that the following inequalities 
hold true 

Ci||ff/ 3 |li 4 exp{-f 2 / 2 } <gf,{y) < Czexpl-fy^} _ (3 22) 

Proof. The results are proven by means of standard super- and sub-solution techniques. We spell 
however the proofs in full details for the sake of clarity. 

To prove ()3.21l) it is somehow more convenient to go back to the variational equation satisfied 
by fp{x) = (\/^£)"i gp{{x - l)/e^), i.e., 

- hfp -hfp + & //3 + ^‘'W{x)fp + ^fl = j,npfp. (3.23) 

The first simple observation is that by positivity of fp and W (x), we get 

-^fp - ilf'p < F (^/3 - fp, 

which, by negativity of the second derivative of fp at any maximum point of fp., immediately 
implies the upper bound 

which in terms of gp becomes, via (13.2011 (here we are assuming that (3 = 


\\gp\\^^=0{l). (3.24) 

V 

In order to prove (13.211) we will provide an explicit supersolution to the equation (13.231) . Notice 
that the first two terms of the equation form the two-dimensional Laplacian, i.e., for any radial 
function /, —A/ = —^dr{rf'). We will use this fact to construct a supersolution in dimension 
two. Let then a > 0 be a parameter independent of e that is going to be chosen later and consider 
the two-dimensional region 

A := Ri_ae 2 ( 0 ) n .4,, = {x G I 1 — < x < 1 — ae^} . (3.25) 


Inside A one has the lower bound 

W{x) > {x —1)^ + 0 ^|a; — 1|^^ > 


2 


+ O > Coa^e^ 
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with Co > 0, so that (13.231) and (I3.2()ll yield 

-iA//3 < - ^^W{x)fp < ^ (Cl - nlCoa^) fp 

where also Ci > 0. If now we pick we get that fp is a subsolution of the following 

differential problem 

-^Af+^Coa^n^e'^f = 0. (3.26) 

To get rid of the inner boundary we now extend fp to the whole ball Bi-ae^ (0) in & smooth (in 
fact at least C^) way. We denote by / such a new function and we require that f{x) = 0 for 
a: < 1 — 2 £^rj and 

-iA/+iCoa2flV/<0, (3.27) 

for any x € Ri_oe2(0). We omit the explicit details of such a construction for the sake of brevity. 
A supersolution to the same problem can be constructed by taking 

/sup(x) :=C'a||/^||ooe-'^('-"') 

with Ca a constant to be suitably chosen: 

-iA/sup + \CoafOf£'^fsup = \ + Co^oa'^Cj /sup > 0, 

if we choose The constant Ca is then used to guarantees that /sup satisfies the proper 

boundary conditions. In order to apply the maximum principle (see, e.g., [El § 6.4.1, Theorem 2]), 
we need that f{x) < fsup{x) on dBi_ae'^{Q), which holds true if Ca > 

/suplag^_^^2(0) = Ga||//3||ooe > //3laBj_^^2(0) ' 


Hence we conclude that / < /sup in the whole Ri_ae2(0), and therefore, using the monotonicity 
of /sup, fp < /sup in the whole region Ri(0) H Ar^. Going back to gp and using (I3.24p . we obtain 
(13.211) in Bi{Q)r\Arf To extend the result to the complementary region, one can use a very similar 
argument with the trivial change — 1 —^ 1 — in the supersolution. 

For the refined estimates (|3.22L we consider the variational equation (13.151) for a < \y\ < g, 
with a > 0 such that and e small enough, which imply 

Up{y) + £^y^v{y) - up = ^y^ - yp + O > ^y^ 


and therefore in that region gp is a subsolution of the equation 

- y - 2(i+e^y) g' + ^ = 0- (3-28) 

As before we extend gp to the whole region \y\ > a in a way and preserving the differential 
inequality satisfied in a < |j/| < g, i.e.. 


-q" 

2g 



Again we skip the details for brevity. 

Now for some C > 0 to be fixed later the following function 

Osupiy) ■■= 


is a supersolution to (13.281) : for e small enough 


--q" 


sup 2(l+e^i/)gsup + °S g^gsup “ ( 4 + 4 (l+e^y) ) 






















CoRREGGi, Dimonte - Third Critical Speed for Rotating BECs 


18 


Choosing the C > ||5/3||ooe°2 to ensure that g^(iLa) < gsup(a), we get the upper estimate. 
Analogously we can choose C > 0 in such a way that 

gsnb ■■= Ce-^y" 

is a subsolution to (13.1511 : first one notes that 

f < /i/3 - +0{e^r]) , 

which follows from the fact that the harmonic oscillator on the real line is bounded from below by 
a/2; then using this inequality in (13.151) . we obtain 


- hdsuh - 2(l+e"j/) gsub + Uf3iy)gsuh + V{y)g^nh + “ Wg^ub 

f - /i/3 + ^gsuh + Ufi{y) - ^y^ + 2 {ilf^y) + 


^sub 


< 


- WgpWlA +0{e^g^) 


9 sub ^ 0: 


if we pick C < Wg/slH 4 - To conclude we use the fact that ^sub goes to 0 as \y\ goes to infinity: 

v 

indeed it is sufficient to observe that there certainly exists a point y > 0 such that gsuhi^y) = 
min{g/3(i7),g/3(-?7)} and 

f g/3(g) |g| < 1 , 

g{y) ■= i g/3(i?) v<y<y, 

[ gpi-v) -y<y< -v, 

is a supersolution to (I3.15L satisfying 5(±j/) > (/sub(±y)- Hence (/sub < g for any \y\ < y, which 
implies the lower estimate (I3.22|) for |y| < 77. □ 


We conclude this Section by stating analogous pointwise estimate for the limiting profile g^\ 

Proposition 3.6 (Pointwise estimates of (/gv)- 
There exists a finite constant C > 0 such that 

Ilggv|l4exp{-fy^} < 5gv(2/) < Cexpl-fy^} _ (3 29 ) 

Proof. The estimate can be proven exactly as (13.221) in ProDOsition l3.5l and we skip the details. □ 


3.3 Optimal Giant Vortex Phase and Profile 

In this Section we investigate the minimization of w.r.t. /3 G R. The main result is the 
following 

Proposition 3.7 (Optimal phase). 

For e small enough there exists a unique minimizer /3* £ R such that 

Ef := inf Ef = Ef. (3.30) 

* /3gR P P* 

Such an optimal phase is explicitly given by 

= + (3.3i) 

where we set g* := gp^ and 

ry 2 rv 1 rv 
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Proof. The existence of a minimizer /3* is guaranteed from the fact that 

Upiy)> 

which implies that = +oo (recall that s > 2). By the same lower bound on the 

potential together with the trivial bound Ef < = 0(1), we also deduce that fii, = O 

In order to find the explicit expression of we first observe that by standard arguments E^ 
is a smooth function of /3 and therefore by the Feynman-Hellmann principled 

dpEf = {gp \dpUp\gp)^ = (^gp (-2^02/ - gp ) ■ (3.33) 


Since /3* is a minimizer, we must have dpE^ 


= 0, i.e., 


£ P* (g* 


ii+Py)'-‘ 


g*) -2flo(5* 

V 


V 


(1+e^ 


g^,) -floe {g* 


V 


(l+e^ 


g*) =0. 

V 


We compute the first and last terms of the expression above: 


g* 


JT+IW 


g.)^ = l + 0{e^), 


(l+e^v) 




g*') — — + O (e^) . 

/ 7] ^ ' 


Indeed thanks to the exponential decay proven in (13.211) . one can easily realize that 

rn 

'-V 

In fact an analogous estimate holds true if g* is replaced with g'^, in particular 

rv 


[ dy \y\'^ gl = 0(1), for any k <oo. 

J-n 


r dyy (g'J =0(1). 

J—n 


(3.34) 


(3.35) 


(3.36) 


(3.37) 


To see this it suffices to integrate by parts and use the variational equation (I3.15|) to go back to 
an expression involving only g^, and there one can use the above estimate. We omit the details for 
the sake of brevity. Notice that at this stage we are implicitly exploiting the bound /3* = 0(e“^), 
which is among the hypothesis of Proposition 13.51 Next we integrate by parts the second term in 
(13.341) to get 


g* 


(l+e^y) 




3*)^ = J_ T^yy* 


rv 

f-r] 


2{i+e‘^y) y^ al 


J -7] 


dy y^ g^gi + ^{l + 0 (e^)) 

I ^ 

where the boundary terms (first term on the r.h.s. of the expression above) can be included in the 


r 

f-T] 


^The notation (■1-)^^ stands for the scalar product in L^([—77,77], (1 + €‘^y)dy). 
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remainder O(e^) if we choose rjo > 2 (see again (I3.21ll i. For the rest we can compute 


/ V 2 n 

dy T+^ 2/^ = --^2 dy (1 + e^y) (y)y*y* 

-T) ^ J —7] 


n ^ 2noe'^P.y + , 

dy-r—757- 9*9* 


'-V 


1 + e^y 


dy (1 + e^y) dy -i [g'^f + ^gt - g-^^gl 


I-7] 


r2 




dy (1 + e^y) y^v{y)dygl 




'-T) 

2 QQy + 


a 


'-1 




dy ■ 


r2 r'^ 


J dy |-i {gif - ^gt + 


-r, l+£2y 

e^(s + 1)17 2£^r2o/3* 


= 


+ 


+ G(e4/3*) + G(e4) = 


— — 2 ~ Q + M* ~ ('5 + 1)^ + 2r2o/3* + O {e^Pfj + O i 


where we have made use repeatedly of (13.3611 and exploited the identity 

{y^v{y))' = da2(s + l)y2 + G(e2|y|3). 

We have also set 




dy (y*) ■ 


(3.38) 


Hence 


9* 


(l+e^y) 




gfj — [— T — sl7 — Q + + 2Ho/3* + O {e^fif) + O (s^)] 


(3.39) 


and plugging this together with (13.351) into (13.341) . we obtain 

J^/3.(l + 0{e^)) + ^[{s- 2)V - Q] + 0(6^) = 0, 

since /r* = T + H + 2Q + 0{e^) (see (13.361) and (13.371) 1. The expression (I3.31|l is then recovered. □ 
Along the proof we have also proven in (13.341) that 


J_ dy (y + ie2y2 _ yj ^ q, 


(3.40) 


which, thanks to the result about /?*, also implies that 


{g*\y\g*)y = 0{e^), (3.41) 

i.e., the profile y* is almost symmetric w.r.t. the origin. 

In fact this latter information can be deduced also by looking at the relation between the 
functional and its minimization and the limiting model From now on we fix /3 equal to 
the optimal value /3*. 

Before discussing this question further we have however to state an useful estimate on y*. 

Lemma 3.1. 

There exists a finite constant C > 0 such that for any y G [—y, rj] 

|y*(y)| < Cy^gfy). 


(3.42) 
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Proof. It suffices to integrate the variational equation (I3.15|) between y > and rj or —rj and 
y < yp^ (recall that yp^ stands for the unique maximum point of 5*): let us assume that y > yp, 
then by positivity of <?* 

^ \9i{y)\ = -i5*(2/) = / {- 2{i+eH) 9* + {^91 + UpAt) + 

pTl 

dt {- 2 (l+e^t )2 + ^91 + UpAt) + 9* + 2{i+e‘^y) 9*iy)- (3-43) 

Now given that the quantity between brackets can be easily bounded from above by Crj^, it only 
remains to use the monotonicity of 5* to conclude the proof. □ 

We are now in position to prove the first result about the energy difference EA — E^'^. As a matter 
of fact this will involve a corresponding statement about the closeness of to in . We recall 
the expressions of the energy functionals 

^r[5] = J dy (1 + e^y) {g’f + UpAy)9^ + e^y‘^v{y)g^ + , 

/ dy{\{g'f + ^y^g^ + ^g^]. 


Proposition 3,8 (Estimate of Ef^ — 

As s —y 0 

EA = E^'' +0{eW), \\9^.-9l4l2=0{eW)- (3-44) 

Proof. We test the two functionals and on suitable test functions. Let us first regularize 
y* outside [—r/,g] to make it an admissible test function for £^: we define 


ytrial (y) 1— Ce * 


9p{y), 

ri (y), 
£2 (y), 


|y| < 9, 
V<y<2g, 
-2y <y<r], 
\y\ > 2y, 


with ri,2 positive smooth functions chosen in such a way that ytriai is at least C^. We also assume 
that both functions ri 2 are also monotonically decreasing. The normalization constant c^, which 
ensures that || ytriai |li2(R) = 1, can be easily estimated: assuming that yo > 2, we have 

Ce = l + OieA, (3.45) 


Since, e.g., 


r2r? 


dy rfiy) < gglAl) = 


Notice also that we need to use (13.4111 to reconstruct the norm of y*: 


'-v 

Now we estimate 


/ dygli^i = cl f dy{l+£^y)gl + cle^ f dy y g^ = + 0{eA- 

J —r) J —71 J —71 


E^ < £^A9t£..y] = cl f dy [\{g'pf + ^y^g} + ^y^} + ^(e"). 
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Thanks to (13.4111 we can easily estimate the error we make by replacing with Up + £^y^v{y)-. 
denoting for short (/) := (5* |/| y*), we have 


(5* -Up{y) - e^y^v{y) 




= { 9* 


a=£=y^(2+£=y) 

2(l+£2y)^ 


+ (2IIoy + - £^y^v{y) g}j = 0 (e^ + ^ 


so that 

i;gv ^Ef' + 0 (e^ + e:^ [y^)) , (3.46) 

where we have also used (13.201) . which in turn requires /3* = G(e“^). 

The trial state for the functional is simply the truncation of g^y, i.e., c^ggv, where now the 
normalization factor can be estimated in this case as 


Ce = l + G(e“), (3.47) 

since by symmetry of ygv 

/ r] poo 

dy (1+e^y) = 1 - 2 / dy 5^^ = 1 + G(e'^), 

-T) J T) 

where we have used the pointwise estimate (13.2911 on g^y and the fact that the integral of a 
gaussian, i.e., the error function, is bounded by the value of the gaussian at the boundary, i.e., 
exp{—c| logep} = 0(e°°) (see, e.g., |AS1 Eq. (5.1.19)]). Then we have 

Ef < Splhggv] = (1 + 0 {e°°)) j dy (1 + e^y) {^(ygv)^ + Up^{y)gl^ + £‘^y^v{y)gl, + ^5gv} 

J- 7 } 

= j dy |i(5;^)2 + +G(e4) = ^;gv + 0(e4)^ (3 43) 

Putting together (13.461) with (I3.48p . we obtain 

Ef'= + O + £^ {y^)) . (3.49) 

Now we decouple the energy E^^: first we bound from below E^ as 

E^^ > J dy (1 + e\) {^(ygv)^ + + ^Sgv} > 

where we have just dropped some positive quantities and used the symmetry of ygv. Then we set 
ygv = uy* for some unknown smooth function u (recall that g^ never vanishes in [—y, yj) and using 
the variational equation for y* as well as Neumann boundary conditions, we obtain 

ES'' >£;f+ J dy {I + e^y) gl + [^y'^ - Up,{y) - e'^y^v{y)j + ^gl{l - 

- f dy M^y*y* + G(e°°). 
J-r] 

Then we estimate 


/ dy (1 + e^y) f ^y^ - Up^ (y) - e 
J-ri ^ 


“^y^v 


2 2 
9*u 


< 


r 

f-T) 


Ce^ f dy|y|y2|l-M^| 

J-ri 

dy (1 + e^y) (^y^ - Up^ (y) - £^y^v{y)'j gl 

<CeW>-\\9ld-u^)\\,.+0{e^f) + e^)^ 
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and by (|3.42|) (notice that the factor 1 + e^y is uniformly bounded from above and below by a 
constant) 


dy u^g*gi 


f dy \l-u^\\g^\\g'J+0{e'^) <Cr]^ f dy \l - u^\gl + 0{e'^) 

J —T) J —T) 

< ||ffJ(l-w^)|L2 +C>(e^), 


which imply by dropping the kinetic term 

>Ef + \ ll5*M'|li2 + i + G (e^ (y^) + e^) 

+ + G (e^ (y^) + eV) • (3.50) 

If we compare what we have obtained with (13.491) . we conclude that 

Ik* -Sgvlll^ =115^(1-w")||l2 (y3)+£V), 


but on the other hand 
rv 


dy y^ g"^ 


f-T) 


r dyy^ {gl-gl^) 

J-T) 




2 _ 2 II 

9 -k ^gV II JJl 5 


(3.51) 


(3.52) 


so that finally 


5J-5gv||r2 =G(eV). 


This proves the second inequality in (13.441) but the first is obtained by replacing the above estimate 
into (13.491) . □ 

It is interesting to remark that a by-product of the proof of Propositionis that (see (13.521) 1 


(5* k^|5*) = 


(3.53) 


which in combination with (13.411) is a strong indication of y* being symmetric w.r.t. the origin 
with a very high precision. In fact this is also made apparent by the estimate of the difference 
5* ~ Hgv. 

The L^-statement in (13.441) can indeed be improved to an L°°-one, showing that y* and ygv 
are pointwise close. The price to pay to have a result in a stronger norm is the restriction of the 
region under consideration to the annulus A,j C Ar^ defined as 


|y e M I ygv(y) > 


(3.54) 


for some 1 / > 0 independent of e. Note that thanks to the pointwise estimates (13.291) and the 
monotonicity of ygv, 

A,J = [-yn,yr,], with » 1. (3.55) 


Proposition 3.9 (Pointwise estimate of y* — ygv)- 
As e —>■ 0 and for any v > 0 


II5* - 5gvlkoo(^^) = G . 


(3.56) 
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Proof. Going back to (I3.5()ll and retaining the kinetic term, we see that we obtain via (13.4411 the 
upper bound 

Wg.u'Wl, = 0{eW), (3.57) 

where we recall that u = ggv/g*- Now let us introduce the set 

€ R I g4y) > , 

so that above inequality together with p.44l) and the bound |1 — m| < |1 — imply 

\w\\%a,) = 111 - = o{eW^^n- 

Then it suffices to use Sobolev inequality in one-dimension: 

||1 - w|Iloo(^,,) < C' (||m'|1^2(_4^) + I|1 - ^IIl2(_4^)) = O . (3.58) 

Finally to obtain the result it remains to observe that Arj C Ari, because in the region where 
g* > 1/(27?'^), by the pointwise estimate, g^v is larger than (1 -|- o(l))/(277'^), which is obviously 
satisfied if gg^ >1/77'". □ 

The above bound shows that inside Arj one can estimate the distance of g* from a perfectly even 
function: for any y € Ar^ 

9*{-y) = 9*iy) + O , 

which is perfectly compatible with the estimates (I3.4ip and (13.531) . 

Another useful consequence of the above pointwise statement is the following 

Corollary 3.1 (Maximum point of g*). 

Let yp^ be the unique maximum point of gi,, then as e —>■ 0 

yp^=0(^e^g^y (3.59) 

Proof The result is a straightforward consequence of the pointwise estimate (I3.5HI1 and the prop¬ 
erties of 5gv (see Proposition [321) ■ 


3.4 Critical Velocity and Positivity of the Cost Function 

From now we fix the phase to be optimal one, i.e., /3 = /3*. The potential function is defined as 


[ dt (1-ke^t) 9,3C;3(i)l/3=/3«ff^ 

^ J —T) 




The main object under investigation is the cost function 

K{y) ■■= yUy) + Fiy), 


(3.61) 


and our main goal in this Section is to prove that it is positive in the bulk of the condensate when 
flo ^ He- To this purpose we will clearly have to investigate the equation (Em and prove at least 
that there exists a positive solution to it. Notice the equation (EJl) involves only quantities relative 
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to the limiting functional fsv jg independent of e. We thus introduce the analogue of (13.6111 
for the limiting case, i.e., the function (j2.19|) 

= kgv + 

where is defined in (12.201) : 


poo 

F^''{y) = -2nJ dttgl^it). 
J V 


We will start by studying the positivity of (12.191) and show that the condition flo > where 
the latter is defined as the biggest solution to (1^ . is sufficient to deduce that K^''{y) > 0 for any 
y GM.. In the second part of the Section we will turn our attention to the cost function (13.6111 and 
prove that the same condition on D,q guarantees positivity of K as well. 

We first observe that F®'' is a negative function vanishing at ±oo: at y = +oo it is obvious, at 
—oo it is a consequence of parity of ggv By this property one can rewrite 

F^'^iy) = 2^0 f dttgl^{t). 

J —OO 

In fact there is another explicit expression of F®'', which can be obtained by using the variational 
equation (13.1711 : 

nOO nOO 

F^^iy) = -^o dtdt{t'^)gl^ = noy'^gl^{y) + 2no dtfg^^g'^^ 

= ^oy'^glviy) h —r / ff'gv [Is'gv - iffiv + 

a Jy 

= -fM7+2T (5gv(y))" + [^oy^ + ^Qois+2) 9l.iy) - (3.62) 

where we have used the exponential decay at oo of g^v to cancel the missing boundary terms. 
Consequently we can rewrite F®'' as 


ia'sviy))^ + + ^oy^ + ^a,ls+2) 9lM - nlfe)] alM- 

The main result about F®'' is the following 

Proposition 3.10 (Positivity of F®''). 

Let rio > 0, then 


(3.63) 


K^^{y) > 0 for any y gW <1 =^ [m®'' “ ^5gv(0)] ■ (3.64) 

Moreover if the strict inequality is verified on the r.h.s., K^''[y) > 0 for any y finite. 

Remark 3.1 (Comparison with |CPRY3) b 

Despite the use of two different potential functions Fi and F2 in JCPRYSf . one should realize that 
\CPRY3l Lemma 3.3] yields the poitwise positivity of a cost function, which is the analogue of K^'' 
in the asymptotic regime flo ^ 1- In fact it can be easily seen that the cost function in \CPRY3f 
is bounded from below by F®'' and therefore the positivity of the latter implies the positivity of the 
first. Hence any threshold CIq one might deduce there must be larger than by definition. 

Proof. One side of the statement, i.e., the fact that the condition F®''(0) > 0 is necessary for the 
positivity of F®'' everywhere, is obviously trivial, so we focus on the other side of the implication, 
namely that F®''(0) > 0 is also sufficient. 
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The core of the proof is to show that 


> 0 for any y G 


> 0 . 


(3.65) 


Indeed if we assume that this double implication is true, a straightforward computation yields 




1 


1 


glio) - 




7rno(s + 2)^s'^' ' no(s + 2) 


5gv(0): 


(3.66) 


since g^v is symmetric w.r.t. the origin and has a maximum at j/ = 0 (see Proposition 13.21) . The 
result is then a trivial consequence of strict positivity of ggv(O). 

In order to prove (I3.65L we first observe that i4rs''(±oo) = 0 and is smooth, so, if there was 
a point j/o where becomes negative, it must be |t/o| < +oo. Moreover as 5 gv, K^'' is symmetric 
w.r.t. to the origin, so it suffices to consider y G R+. The derivative of is easily computed 
from the expression (12.1911 : 

= Ogvg'gy + 2f2o2/ffgv, (3-67) 

and one immediately has that = 0, i.e., has a critical point there. Whether it is a 

minimum or a maximum depends on s and ilo; but as we are going to see this does not matter. 
We can in any case compute easily the second derivative of exploiting once more p.l7l) : 


= (ffgv)^ + 4noyffgv5gv + 2 [^a'^y'^ + flo - + ^5gv] ffgv 


(3.68) 


Then we prove the crucial property of /•f®'': suppose that has a maximum at j/i > 0 and 
then a minimum at 2/2 > yi, then 


K^^{y2) K^'^jyi) 
glAy^) ~ gt,{yi) ’ 


(3.69) 


and in particular K^'^{y 2 ) > 0 if > 0. 

To conclude the argument once (13.691) is proven, it is sufficient to observe that iG®'' has a critical 
point in ?/ = 0, which by parity must be either a maximum or a minimum: if it is a maximum, 
then (|3.69p shows that at any minimum point 2/2 > 0, i^®'"( 2 / 2 ) > 0. Notice that it does not matter 
whether if®'' has a single or multiple minima, because any minimum after the first requires the 
presence of a preceding maximum point, where is larger than its first minimum and therefore 
positive. If on the opposite K^'' has a minimum at the origin, then it means that there must be a 
maximum at some 2/1 > 0, where obviously > K^''(0) > 0 and we can repeat the argument 

for any minimum after 2 / 1 . 

Let us now prove (13.691) : we assume again that it'®'' has a maximum in 2/1 > 0 and a minimum 
in 2/2 > 2/1- Then it must be K^'{yi^ 2 ) = 0, i.e.. 


ggv{yi,2) — -2fio2/l.2ffgv(2/l.2)- 

Moreover replacing this condition in (12.191) and (13.6811 , we get 

K^{yi,2) = klo^yl2 + 1- n7(f+2) (m®" - ^5gv(2/l.2))] 5gv(?/l.2), 


Moreover 
which implies 


A"®'"'(2/i, 2) = [fIo('S - 2)yi,2 + 2fIo - 22^®'' + f 5gv(yi.2)] 5gv(yi,2)- 
i^®""(2/i) < 0 < if8""(2/2), 


(3.70) 

(3.71) 

(3.72) 

(3.73) 


- 2 ) 2/1 + 2f2o -2g+ ^gl^iyi) < ng(s - 2)yl + 2^0 - 2/r + f 5gv(2/2) 


(3.74) 
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and using this inequality in the expression of if®''( 2 / 2 ), we obtain 


KS''(y2) 

ffiv(y2) 


> flo 


s-2 


s + 2^1 7rfio(s + 2} 

KS-(yi) 


(3gv(y2) - Sgv(2/i)) 


1 


“1“ 2) 




ffiv(yi) 7rOo(s + 2) ^ 


(ffgv(2/i) - ffgv(2/2)) > . V ^ (3-75) 


ffiv(yi) 


because by hypothesis 2/2 > 2/i ^-^d g^v is decreasing. 

Notice that as a by-product of our analysis we found out that can have no global minima, since 
infygRif®''(2/) = 0 and therefore at any such minimum point 2 /O 1 we would have if®''(2/o) = 0, but 
this clearly contradicts (13.691) . Hence if the inequality on r.h.s. of (13.641) is strict then if®''( 2 /) > 0 
for any finite y. □ 

Proposition 13.101 introduces the equation (12.71) . The next step is obviously to prove that such an 
equation as at least one solution; 

Proposition 3.11 (Equation (12.71) 1. 

The equation dZIl) 


Hn — 


s 2 


has at least one solution Hq > 0. 
Proof. Let us first set 


G(Ho) := 


s 2 






27T 


so that (EZl) reads Hq = G(Ho). We will show that G(Ho) is asymptotically smaller than Hq (resp. 
larger) Hq for large (resp. small) Hq. 

Let us first consider Hq ^ 1: using the trivial bound ||(?gv|l 4 < llUgvIl^ and the definition of 
2 ^®'', we get 

G(Ho) < 


If now we plug into f®'' as a trial state the ground state of the harmonic oscillator ft-osc = + 

we easily obtain 

E®'' < ^noyTT2 (1 + , 

so that 

G(L!o) < (1 + G( 0 -^/ 2 )) < Oo, 

if Ho » 1, because s > 2. 

On the other hand for small Hq, thanks to the estimate (13.181) and again the definition of /r®'', 
we have 

G(Ho) > ^E®''. 

To bound from below E®'' for small Hq, we can simply drop the kinetic term to get 

E^'' > inf f dy ^ 

llp|li=i Jr 

i.e., a TF-like functional. By scaling we immediately obtain 


as Ho —?► 0, so that 
for Ho small enough. 


inf . / dy {iaVy+ = GHof^ 

p|li=i Jr 

G(Ho) > GHq/^ > Ho, 


□ 
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In the above Proposition we have not investigated the uniqueness of the solution. We indeed 
expect that such a solution is in fact unique, but without a proof of this fact, we have to choose 
Oc equal to the largest possible solution. 

Now we turn our attention back to the cost function K: as K^'' it is given by the sum of a 
the positive density and the negative potential function F (see next Proposition I3.12() . In 
addition g* is monotonically decreasing for y > = o(l) and, like F, almost symmetric. Close 

the boundary of the interval [—77,77], <7* gets extremely small (in fact exponentially small in e) but 
F vanishes identically at ±77. In conclusion it is clear that the overall positivity of K should then 
emerge from a very delicate balance between the two opposite contributions. 

We first state some simple properties of F collected in the following 

Proposition 3.12 (Properties of F). 

The potential function defined in (13.601) is such that 

F{y)< 0 , for any y £[-r],r]], (3.76) 


F(±77) = 0. (3.77) 

Proof. One of the identities (13.771) is trivial, the other is a direct consequence of (I3.40|) . In order 
to show that F is negative everywhere we compute the derivative 

F'{y) = T+^ (2^02/ + - I3*£^) gliv): (3.78) 

and it is easy to verify that because of the first term F'{—rj) < 0 while F'{r]) > 0. Moreover 
F' vanishes at a single point yp = G(e^), where F has a global minimum. Hence it is negative 
everywhere in [—77,77]. □ 

A very crucial piece of information about the potential function formulated in the next Propo¬ 
sition is an alternative expression of it, which relies on the variational equation (13.151) and is the 
analogue of (I3.62|) for F®'". 

Proposition 3.13 (Alternative expressions of F). 

For any y G [—77,77] the potential function F admits the following alternative expressions 


F{y) 

where 


^0 , ,.2 
- t[9J + 


2Hn 


1 2 2 , ^ 2 

y + - y* 


R±iy) = 0 {e^rf)gl{y), R± = 


^2 \R+iy) + R+^ ify>yy., 
\R-{y) + R-, ify<yy^, 

(3.79) 

-^2oy^5*(±y) (1 + 0(1)). 

(3.80) 


Proof. We consider only the case y > yp, since the other one is analogous. The key ingredient of 
the proof is an integration by parts, exactly as for (I3.62|) . We spell all the details nevertheless for 
the sake of clarity. Thanks to the vanishing of F at 77, we have 


F{y) = f dt (-2flot - -b/3*e^) 5* = / dt gf dt {-^ot'^ - 

Jy Jy 

(-H0772 - gf{y) - - \^o£^y^ + /3*e^y) gl{y)+ 




+ dt [—flot^ — 




= -noy^gl{y){l + o(l)) + + 0 {e^y'^)) gl{y) 


- 2 


/ g^^gi. (3.81) 

J V 
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We now rewrite the last term by reconstructing the potential (7 and using the variational equation 
(13.151) : since 

- = ^(1 + eH)Up^ (t) 

+ - 2)/3*t + , 

we obtain 

/■’* 1 40n n 

“ ^ y {l + e'^y) “ \VLoe^t^ + I3^£^t) g^g'^ = — J dt (1 + e^t) Ug^ {t)g*gi 

+ 2 £^ f dt (^CeH) “ iI 2 o(s - 2)/3*t + 3 * 5 * 

= ^ f di (1 + eH) UpSi)g*g* + 0{£^g^)gl{y), (3.82) 

« Jy 

where we have used the bound (13.421) and the monotonicity of (/* for y > yp^. The first term on 
the r.h.s. can be rewritten by means of (13.151) : 


dfin r 


2 Q r 

dt (1 + eh) Up^ = -Zr dt (1 + eh) i (gi)^ - ^gt + y^g^ 

- a ^ Jy 

+ —^ / dt ( 5 ^)-^ / dt (1 + e^t) t^v{t)g^g^ = 

= ^ (1 + [-^dtiv) + f^*9Uv)] - ^ (1 + e^y) 5 (fl*(y))^ - ^y*(y) + ^^*9l{y) 


2 Hq£^ 


dt 


1 / \2 I 1 4 g 

2 (y*) + ^y* - y*y^ 


4ftoe^ 


dt (1 + £^t) t^v{t)g^,gl 


= ^y*y*(y) - ^ 5 (yl(y))^ - ^y^y) + M*y?(y) + c>(£^y^)y*(y)- 

Putting together the above estimate with (|3.8ip and (|3.82p . we obtain the result. □ 

Thanks to Proposition 13.131 the cost function K can also be expressed as 


h 2 / X 2 f 2 o 


K{y) 


r^o . / .2 , 2 rio 

2 \9ir) 4“ 2 


a 


2..2 


4fto + 2 “ y 




gl 


R+{y) + R+, 

R-[y) + 


if9>w., (383, 

if y < y/ 3 . • 


This alternative expression will play an important role in the proof of its positivity, exactly as for 
. Another important ingredient of the proof is also the closeness of K to as £ —>■ 0: 

Lemma 3.2. 

For any flo > 0 and y € Arj 

K{y)-K^hy) = 0{e^\\oger). (3.84) 

Proof. The result is a direct consequence of the pointwise estimate (13.561) . □ 

The above Lemma in combination with Proposition 13.101 might seem to give also the positivity 
of K inside A,,. However this is not the case because, although we proved that AT®'' is positive 
on the whole real line, we did not provide any lower bound to it. In fact even by just looking at 
its minima, one could conclude from (13.691) that K^'^{jj 2 ) > KS''{yi)g^^{y 2 )/g^^{yi), where 1 / 2 ,yi 
are the positions of the minimum point and the preceding maximum point (consider for simplicity 
the half-line R+). Now even if K^''{yi) > C > 0 as it occurs for instance at the origin, the ratio 






































CoRREGGi, Dimonte - Third Critical Speed for Rotating BECs 


30 


between the densities can become extremely small in Ar). In addition to that the inequality holds 
true only for the minima of and it might be that it has no minimum inside Ar; or Arj , in which 
case we only know that it is positive there, but without any meaningful lower bound. 

In fact we will be able to prove positivity of K only in domain strictly smaller than Ar/^ because 
of the additional constant terms R± in (I3.83|l . irrespective of their smallness. We thus set 

•A> ■= {y e [-r],v] I gliy) > uiax{gl{r]), gli-g)}} . (3.85) 

By monotonicity of g* for large y is easy to see that A> = [—y_,t/+] with ?/± —>• oo as e —>■ 0. 
Notice also that 5 * is very small at the boundary of A>, although not as small as g*(? 7 ). 

We can now state the main result of this Section: 

Proposition 3.14 (Positivity of K). 

If^o > us £ -^ 0 ; 

Kiy)>0, foranyy€A>. (3.86) 

Proof. As in the proof of positivity of AT®'' in Proposition l3.10l the key idea is to show that positivity 
at the origin is indeed sufficient to get the result. This in turn is easily inherited from positivity 
of whenever fig > fic, via (13.841) : by (13.661) 

a:®''(o) >C>0, 

but 

A:( 0 ) - A:«''( 0 ) = ©(e^l loge|“), 

and thus Ar(0) > C > 0 for a possibly different constant C. 

The rest of the proof follows the same line of reasoning of the proof of Pronosition ld.lOl There 
are however two complications: first we have two alternative expressions of AT in {yg.^,g\ and 
[—g,yp^] respectively. Recall that yp^ = o(l) denotes the unique maximum point of g*. Second 
the presence of the constant terms R± in (13.831) is very annoying and in fact it is responsible of 
the restriction to A>. 

In order to handle the first issue it is sufficient to take into account the two intervals [yp^ , g] 
and [—g^yp^,] separately and use a different expressions for K (see (13.831) '). 

The second issue on the other hand leads to the introduction of the modified cost function 

K{y):=K{y)-6,gliy)-R± (3.87) 

for some 

0 < 4 < < 1 (3.88) 

to be chosen later. Here we have used a compact notation to mean that we subtract i?+ (resp. 
R-) in [yp,,g] (resp. [-g,ypj). 

Now we observe that if Hq > He 

KiypJ = K{0) + o(l) = KiO) + o(l) > 0, (3.89) 

thanks to the pointwise estimate (13.841) and since ||g*||ge) ^ G. It is interesting to remark that this 
is the only point in the proof where we use the condition Hg > He, although several later estimates 
are affected by this one. Moreover at the boundary of the domain we have 

K{±g) = {^-S, + 0{e^g^)) gl{±g) > 0. (3.90) 

Therefore in order to exclude that K becomes negative, it suffices to prove that it is positive at 
any possible global minimum point —g < y^ < g. 
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We claim that, for flo > He, any global minimum point of K must satisfy the condition 

|2/m| » 1. (3.91) 

The reason is the pointwise estimate (13.8411 and the observation contained in Proposition 13.101 K 
and are pointwise close and therefore we would have 

K^''{ym) < uAn[K{-ri),K{ri)\+o{l) < \ram.[gl{-ri), gl(jq)] + o(l) = o(l), 

which in turn implies K^{ym) < 0 since is independent of e. For flo > He this contradicts 
the statement of Proposition l3.1()l ^ 

The keyjDoint in the proof is the following property: let ym be a point where K reaches its global 
minimum Kijjm) < (otherwise there would be nothing to prove) and yM any maximum 

point of K such that j/m < 2 /m, if J/m > 2//3*, or 2 /m > 2 /m in the opposite case 2/m < 2//3*- Notice 
that such a maximum needs not to be the global maximum but its existence is a consequence of 
smoothness of K and the inequalities 

.^(± 77 ) < KiyyJ, Kljjm) < Kijjfi,)- 

Then we are going to prove that 


K{y^) K{yM) , , 
glivm) - gliVM) ^ ^ ^ 


(3.92) 


Now suppose that this is true, then we can pick a maximum point yM of K where K{yM) > 
> G > 0. In addition it must be 


2/m = G(1), (3.93) 

because K{y) < Cg'^{y) and the decay estimate p. 22 l) or the pointwise estimate (13.561) imply that 
K{y) = 0 ( 1 ), if I 2 /I ^ 1. Hence by the lower bound (13.221) gi,{yM) > G > 0, (I3.92|) yields 

KiVni) > gliym) (c~^K{yyJ + o(l)) > CogHym) > 0 (3.94) 

for some Gq > 0. In fact we have obtained something more: for any y € Ar; either K{y) > 
mm{K{—r]),K{r])} > 0 or 

Kjy) > Kjym) ^ gliym) „ 
gliy) - gliy) “ gliv) ' 

Either way K is positive everywhere in [— 77 , 77 ]. Moreover the positivity of K implies that 
K{y) > R±+6egl{y) > 0 , if gl{y) > 

for any y G Arj. If now we restrict the inequality to .4> and we choose, e.g., Se = the estimates 
(13.8011 imply that inside 

gliy) >ifTa&y.{gli-ri),gl{g)) Cif xn&-K{gli-r]),glig)) 

so that Kiy) is strictly positive for any y G ^>. In fact a closer look to the chain of inequalities 
reveals that we have proven something more, i.e., for e small enough 

Kiy) > v~^gliy), 


for any y G M>. 


(3.95) 
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We have now to prove (13.9211 . Recall the assumption: we have a global minimum of K at ym 
and a maximum at yM, which is on the left (resp. right) of ym, if ym > y/ 3 * (resp. ym < yp^)- The 
idea is the same used in the proof of Proposition 13.101 the derivative of RT, i.e., 


K'{y) 


(1 - 4)5* + ( 2 floy + 


/3*e^) 5* 


5*, 


(3.96) 


must vanish both at ym and j/m and therefore 

1 


9*{ym,M) — 


(1 - (5e)(l + e^5m.M) 


(2flo5m,M ~t“ Um.M Z^*^ Q'kijJni.M') • (3.97) 


The second derivative of K can be computed as well: 

K''{y) = (1 - 4) {g'S + 4b!oy5*5; + [(1 - 4) f 5*" - 2/i*) + 2flo + 0{e^y^)\ gl (3.98) 


so that at any extreme point of K, one has 

K''{ym,M) = [^o(® ~ ‘^)y‘rn,M + 2f2o + ^gl{ym,M) ~ 2^* + o(l)] 5*(5m,M), 
where we have exploited the condition (I3.88|) . Similarly by (13.831) we get 

K{ym,M) = + 2 ~ 7T!7(s+2)5*(5m,,M) — 5* + o(l) gi,{ym,M), 

and a direct comparison between (13.991) and (13.10011 yields 

^ 2 9i^{ym,M) TC {jJm,M^ 


9l[ym,M) 5 + 2 7rno(s + 2) ^o{s + 2)gl{ym,M) 

Now this is the key identity because by assumption (recall also (13.911) 1 

K"{yM) < 0 < K"{ym), 9*{ym) < 9*{yM), 


+ 0 ( 1 ). 


(3.99) 


(3.100) 


(3.101) 


so that _ _ 

K{ym) K{yM) , . 

9l{ym) - 9l[yM) ^ 

i.e., (13.921) is proven. Note that the fact that we have two different explicit expressions of K for 
y > yp^ and y < yp^ did not affect the proof, because the difference between the two expressions 
is o(l) and therefore can be included in the error term. □ 


4 Energy Asymptotics 

We attack in this Section the proof of Theorem l2.21 which will imply the main result of the paper. 
The result is obtained by combining upper (Proposition |TT|) and lower fProposition 14.21) bounds 
on 

4.1 Upper Bound 

The upper bound on is stated in next 

Proposition 4.1 (GP energy upper bound). 

A.S s —)• 0, 
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Proof. The proof is rather simple because it is sufficient to test on suitable trial function of 
the form 

^triai(x) := ^ gtriai ( exp {i [O + /3*J 9} , (4.2) 

vZtts: '' ' 

where Atrial coincides up to a normalization constant with g* within Ar^ and is suitably regularized 
outside. The calculation is rather straightforward and we omit it for the sake of brevity. Note 
that the remainder 0(1) is entirely due to the fact that the phase + /3* might not be an integer 
number. Otherwise one would obtain a much better error term O(e^). □ 

4.2 Lower Bound 

A lower bound for matching the upper bound of Proposition |TT] is formulated in next 

Proposition 4.2 (GP energy lower bound). 

If flo > flc, as e —>• 0, 

pgv 

>^ + 0(e°°). (4.3) 

Proof. We first restrict the integration in the GP energy functional to the domain Arj (recall its 
definition in (ESI): to this purpose we just have to observe that all the three terms in the GP 
energy functional are pointwise positive and thus we can simply drop their integrals outside Ari- 
Of course is not normalized in L^{Ari) but the exponential decay proven in Proposition 13.31 
guarantees that 

by taking tjq large enough. 

The first step in the proof is a splitting of the energy, in order to extract the leading order term 
Ef /e'^. This is now rather standard and we do not spell all the details of the computation. We 
just note that one sets 

?/>^P(x) =: -^u {x, d) g* (^) (4.5) 

y/27Te 

Since ri + /3* needs not to be an integer, u is not single-valued in general, but 

u{x, d + 2fc7r) = d), (4.6) 

for any fc G Z. A part from that u is finite for any x S Ar; thanks to the strict positivity of g*. A 
long but simple computation using the variational equation for gives 

E^^ >^ + ^+0{en, (4.7) 

where the inequality is mainly due to the restriction of the integration domain and setting y = 
1 + e^x for short 


£[u]= j dx g 2 (y) |i |Vu|Va.j„ + _ |y| 2 ) 2 | 

d .jA-Y) 


a{x) := 



5 


(4.8) 

(4.9) 


and the superfluid current is defined in (12.161) . The rest of the proof is devoted to prove that 

£[u] > C>(e°°). (4.10) 
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In order to exploit the cost function trick mentioned in Section 12.11 and the positivity of K 
proven in Proposition 13.141 we need to restrict again the integration domain in £[u] to C A-q, 
where 

Af^ = {x e I |1 - |x|| /e^ e ^>} . 

The only annoying term is the only one which is not positive, i.e., the second one in (14.81) : 


/ dx5r2(y)a-jt 

JAr,\Al° 


< l|a|lL~(^„) 


’Aq\Ai^ 


dxg*(y) |m| |Vi,u| 




T? \'^> 

-4_4 IL/GP 


(4.11) 


where we have used the bound ||V'!/'®^|| raa/ a ^ ^ Ge following from 

II 11 \.”4,t7 ) 


I|V^ 1 L<G 




lloo)> 


(4.12) 


which can be proven from Gagliardo-Nirenberg inequalities exactly as in [CRYl Lemma 5.1]. How¬ 
ever the lower bound (13.221) easily implies that if we set Ay =: [—y-,y+], then 

y± = ? 7(1 + o(l)), 


so that 

(4.13) 

again by (13.191) and the arbitrariness in the choice of r]o. Hence (14.111) yields an error which can 
be made smaller than any power of e and we get the lower bound 

£[u] > _^^^dx gl{y) |i |Vu|Va.j„ ^9l{y){l - \u\‘^y'^ + 0{e°°). (4.14) 

We can now finally integrate by the angular momentum term by using the potential function 
F defined in (13.601) : it is trivial to verify that 

2^2 (^) = -d^F (^) (4.15) 

so that 



dx5*(2/) a- 


1 /• 27 r y+ 

- d-d dx dxF ^[iu{x,'d)dsu*{x,'d)] 

2 Jo Jl-e^y- 


1 

2 


fl+e y+ 


dd 




dx F 5? \idxu{x,'d)d.yu*{x.-d) +iu[x,'d)d1,gu*{x,'d)\ 


1 

2 


r2-n 


d-d 


F 


^ [iu{x, 'd)dsu*{x, -d)] 


i+Ay+ 

l-e^y- 


(4.16) 


The boundary term can be easily proven to provide an exponentially small correction: consider, 
e.g., the term at 1 -I- e^y+, since |F( 2 /±)| < C-q^g^F-q), one can reconstruct a term, which can be 
bounded exactly as (14.111) . The result is an error 0{e°°). The rest is integrated by parts once more 
but this time w.r.t. -d: 


1 /.27I. Al+e^y+ 

- d-d dx F {^)^[iu{x,-d)df^u*{x,-d)] 

^ Jo Jl-e^y- 


1 ^27r |•l+E‘^y+ 

- dd dx F (^^^^)^[id^u{x,d)dxu*{x,d)], 

2 Jq Jl-e^y- 
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where the vanishing of boundary terms is due to the periodicity of u*dxU and its complex conjugate 
(compare with (14.61) 1. Altogether we have thus obtained that 

[ dxgl{y)a-iu= f dx F {^)^[iWxu{x,'d)W^u*{x,d)] + 0{e°°) 


> - 


lAl^ 

f dx\F(2^)\\yuf + 0(e°°) = -f dxF(2^)\Vuf+ 0(e°°), (4.17) 

Jai° Jai° 


and therefore 

£[u] > dx (^) |Vu|^ + ^5?(2/)(l - + 0[en 

[ dxgl{y)\Vu\^+Oie°°)>Oie°°), (4.18) 

Jai° 

thanks to Proposition [TH] and in particular (13.9511 . □ 


5 Giant Vortex Transition 


In this Section we prove the results regarding absence of vortices and total vorticity of the conden¬ 
sate. 


Proof of Theorem \2.1\. Combining (14.181) with (14.711 and the upper bound proven in Proposition 
ITTl we get 


/ dxgt{y){l-\u\^f = 0{e% 
JAiP 


(5.1) 


which already means that \u\ can not differ too much from 1. To deduce the pointwise estimate of 
Theorem O we need to combine this with an estimate of ||Vu||j^. 

As in |CRY1 Lemma 4.3] we obtain from (11.1611 and (13.1511 the following variational equation 
for u: 

- - ^g'^dxU - ig^^a ■ Vu + ■:^gl (^\u\^ - u = - ^/i*) (5.2) 


which yields (recall the definition of .Abuik C in (12.811 ) 


IIAull, 


= Mbuik) A C 


^g\\Wu\\ 


(«4.bulk) 




Now using the elliptic estimate (I4.12p we conclude that 

= O (£-“'f‘+*) . (5.3) 

Suppose now that it exists Xq £ -Abuik such that |u(xo) — 1| > loge|*', for some 6 > 0 to 
be chosen later. Then from (ESI) we get that 

\\u\ - 1| > loge|^, for x £ (xq) (~1 A, 

with p = e®/^| loge|^“^“^, and 

0{e^) = f dx gtiy) (l - \u\^f > Ce^\\oge\^^-^‘^-\ 

J ^buiknBe(xo) 

which is a contradiction for all 6 > 4a — 1. □ 
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We now focus on the proof of Theorem [531 and for later purposes we state a useful Lemma, which 
is the analogue of [CPRY31 Lemma 3.5]: 

Lemma 5.1. 

Let fig > fic and R be a radius satisfying R = 1 + O (e^), then 

|deg(u,5Rfl)| = G(l) (5.4) 

Proof. We use a smooth radial cut-off function x with support in [i?, i?] such that x(i?) = 0 and 
x(i?) = 1, for some radius 

R = R — cs^ 

with c > 0. We also require that |x| < 1 and |Vx| = 0{e~‘^). Then by Stokes formula 


deg {u, OBr) = 


'OBr 


_ 1 /" 
u ) tt 


daxiR) 3 


dx V X ■ ^ I - 

u 


Vu 


(5.5) 


Therefore 


\deg{u,dBR)\<^ f ( 5 g) 

£ JbrXBs m ^ “ 


where we used that ||1 — = o(l). Now the result proven in Proposition 13.101 in fact 

says that for any x G Br \ B^ and for flo > He, there exists a constant C > 0 such that 

^gv 


and thanks to (13.841) the same inequality holds true for A', i.e., 


^(^)>G>0, for any X G Rfi \ (5.7) 

Going back to (14.181) this yields 

Oie^)> f dxiG(^)|Vu|2> / dxiG(^)|Vu|2>G||Vu||^.(g^\g ), 

JBR\Bij^ 

which gives the result once plugged into (15.61) . □ 

We are now in position to complete the estimate of the winding number of : 

Proof of Theorem \ 2. S[ We follow [CPRY31 proof of Theorem 1.5]. The positivity of on OBr 
is guaranteed for any radius R = 1 + O (e^) thanks to (12.91) . A simple computation shows that 

deg , OBr) = H -h /3* -h deg (it, OBr) , 


which yields the result in combination with (1^ . 


□ 
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